NPS ARCHIVE 
1966 

RADECKI, R. 



SHALLO'i ’./AT.R T/^ROACOLS'iIC COr. [HICATICRS 



by 



RICKARD AITTHOIW RAD3GKI 



MSSACHUbi^TIS lASTITUTD OF TECHNOLOGY 

Thesis 

R135 

MY, 1966 






DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA9394«101 



1* « Qji'e Schor/f 
« l!# *v<itiA 






SHALLOW WATER HYDROACOUSTIC C OKrlUNICATIONS 



by 

Richard Anthony Radecki 

& 

B.S., United States Naval Academy 

1959 

SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS 
FOR THE DEGREES OF NAVAL ENGINEER AND MASTER 
OF SCISN3E (SLEGTHIGAL ENGINEERING) 

at the 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
June , 1966 



Signature of Author 

Department of Naval Architecture and Marine Engineering, 

May 20 , I966 



Certified by 



Thesis Supervisor 



Accepted by 

Chairman, Departmental Committee on Graduate Students 



» y 







y T-f- -xduate School 




Kfa& 

^Al)PClt(, 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5101 



SHALLQ'i WATER HYDROACOUSTIC COMMUNICATIONS 



by 

RICHARD ANTHONY RADECKI 



Submitted to the Department of Naval Architecture and 
Marine Engineering on I-lay 20, 1966 in partial fulfillment 
of the requirements for the Toaster of Science degree in 
Electrical Engineering and the Professional degree. 

Naval Engineer. 



ABSTRACT 

Hydroacoustic communication has not been used ex- 
tensively at ranges in excess of ten miles. A study of 
the shallow -vfater and surface ducting channels indicate 
that reliable sixty word per minute communication at ranges 
of fifty miles is possible under a variety of conditions, 
with reasonable power levels and equipment complexity. The 
channel can be approximated by a Rayleigh fading signal with 
additive Gaussian noise. Five frequency multiplexed binary 
channels are compared with a fifteen channel system utilizing 
frequency diversity, and with a fifteen channel system util- 
izing coding. ^ method of calculating the diversity error 
for components with unequal signal and or noise levels is 
worked out. Frequencies must be above 900 cps if shipping 
noise is significant, and belovj 15OO cps to minimize at- 
tenuation. '-^'o optimize several transmitter-receiver para- 
meters, such as the signal interval, more long range ex- 
perimental data is required. 
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INTRODUCTION 

Except for a fevr very short range or very lovj rate 
systems, the cominunications capability of ocean acoustics has 
not been exploited. Voice modulation has been used extensively 
for ranges of a few miles and CW at slightly longer ranges. 
Single explosive charge signals have been sent over a thousand 
miles via deep ocean channels. Aside from sonar these few 
examples cover most of the communications development. Vjith 
the increasing interest in nuclear submarines, oceanography, 
and the industrialization of underwater real estate, the use- 
fulness of a reliable system which does not depend on radio 
or cables is apparent. 

In this thesis a single such problem will be con- 
sidered, that is, communication between two stations in rela- 
tively shallow water using transducers near the surface. A 
range of fifty miles and a data rate of sixty words per minute 
will be used to evaluate various systems. If such a system 
were feasible it would provide a good link between ships ope- 
rating in the same geographic area, or betx’reen an oceanographic 
vessel and its remote data collecting stations, or any number 
of other purposes. 

Part I of this thesis will be a simplified description 
of the shallovr water acoustic channel, which will point out 
the parameters important to a communications system. These 
include ambient noise levels, transmission losses, and channel 
stability, •‘■'hese parameters will be used to arrive at a 
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random process model of the channel. 

In Part II this model vjill be used to develop a com- 
munication system. First a single binary channel will be 
developed. This will then be extended as required to obtain 
the desired data rate at a range of fifty miles with an ac- 
ceptable error rate, ^ince there are many vrays to handle the 
problem only one or two representing current technology will 
be discussed. 

In Part III the feasibility of hydroacoustic communica- 
tions under various ocean conditions will be evaluated. The 
assumptions used in the first two parts and their effects will 
also be evaluated at this point. 

Part IV will summarize the results and indicate areas 
requiring additional research and or development. 
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Part I 

SHALLa-i -VJATER HYDROACOUSTICS 

1. General . In shallow vzater there are two principle methods 

of acoustic transmission. I’he first type is the surface channel. 
In this channel a lar/re portion of the acoustic energy is 
trapped in a duct near the surface. Transmission in this 
duct is by surface reflection and upward refraction. The second 
type is transmission by both surface and bottom reflections. 

The particular type of transmission which vrill occur is pri- 
marily dependent on the vertical profile of the velocity of 
sound in the ocean area. The velocity of sound in sea water 
varies primarily with temperature and to a lesser degree with 
pressure and salinity. Over long ranges it is reasonable that 
both types of transmission may occur for part of the distance. 

2. The Surface Channel . The density of viater decreases with 
temperature. This means that in general the water near the 
surface is warmer. However, if a cooling condition exists, 
the surface water may be cooler or at the same temperature 

as deeper water. This inversion causes the surface channel to 
exist. The velocity of sound increases with increasing pres- 
sure and with increasing temperature (1). Therefore, an 
inversion or isothermal layer will have a positive velocity 
gradient, i.e., velocity of sound increases x-xith depth. Ac- 
cording to Snell's law, acoustic waves in this layer will be 
bent upvxard, reflected from the surface and bent upward again. 
Figure 1 is a diagram of this type of channel. It has been 
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Figure 2. Spherical and cylindrical spreading losses. 
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shovm and experimentally verified, that a minimum layer depth 
is required to trap a given frequency signal (2). This depth 
is : 

Iffiin = 2.7 X 10^ (fg)“^^^ 

f = frequency (cps) 

g = duct velocity gradient (sec.~^) 

For a 1000 cps signal, and a 0.020 meter/second/meter velocity 
gradient this minimum layer depth is 30.4 meters. 

There are several types of transmission losses which 
must be considered in this channel. These are: spreading 

losses, reflection losses, and absorption losses. The spread- 
ing loss will be broken into spherical and cylindrical com- 
ponents, The energy leaves the source initially as a spherical 
wave. Beyond a certain point the energy is essentially trapped 
in the layer and the spreading loss becomes cylindrical. For 
purposes of calculation the transition point t-iill be taken as 
half the horizontal range of a ray leaving the source at the 
layer depth horizontally and ;)ust reaching the surface. This 
is shown in figure 2. 

The second type of loss to be considered is the surface 
reflection loss. The reflection coefficient for the viater- 
air interface is approximately minus one. This is true for 
calm water, but waves cause a certain amount of scatter loss. 
This loss coefficient has been determined experimentally (3)» 
The follovxing formula describes the loss function for a fre- 
quency -waveh eight product less than 3: 



/ 
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-10 (1-0. li^9 (fH)3/2) db per bounce 

f = frequency in kc/sec. 

H = waveheight in meters 

'The third type of loss to be considered is absorption. 
This takes into account primarily the thermal and viscous 
forces which are otherwise neglected in deriving the vrave 
equation of acoustics in a liquid medium. For the ranges of 
interest this loss is negligible belovi frequencies of 2 kc. 

A suitable formula for absorption losses which is accurate 
from frequencies of 2-25 is (4) ; 

’^ 2 .= ^ — £_ + 3 *35 X 10~® f^' nepers per meter 

7750 + f2 

this gives a transmission loss of: 

L =7,280 X 103x db/km. 

If it is assumed that the velocity of sound does not 
vary in the channel except linearly vfith depth, some rather 
simple equations can be derived xMhich specify the geometry 
of the path in depth and range. Sef erring to figure 3, 
these formulas are (5 ) J 

d = transducer depth (meters) 

L = layer depth (meters) 

R = horizontal range (meters) 
a = velocity gradient (sec.“l) 

90 ° - angle of incidence 

Y = vertical depth used for calculation (meters) 

Cq= lowest velocity of sound in the layer (meters/sec.) 

( 
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Velocity profile. 



Typical transmission paths. 






Figure 4 . The SOFAR channel 
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Y = CO , 

~ (cos cos -Oi ) 

a cos •©; 



R 



Co 

a c os 



(sin sin ■©■ ) 



cos 

Co c o+aY 

If we now let Y = L (transducer at the surface), and ■©: = 0 

then 2R is th^maximum distance betvj’een reflections, and the 

minimum number of reflections over a distance D is: 

minimum num/oer of reflections = ^ 

2R 

Since the vjater depth does not enter any calculations 
for the surface channel, it is also a suitable deep vrater model 
as long as it is uninterrupted. 



3 . 1*he Bottom Reflection Channel . The bottom^ reflection 
channel only exists in shallow water. In very deep water the 
temperature reaches a minimum at some intermediate depth and Is 
isothermal down to the bottom. This gives rise to a minimum 
velocity of sound, i^hich generally occurs at depths of about 
two thousand feet. Sound waves are reflected upward below 
this depth, and downvjard above this depth, ’-^'his is the SCPA.R 
channel shown in figure 4. It will not be considered in this 
thesis. For the shallow vzater channel, ocean depths up to 
600 meters will be considered. This includes about ten per 
cent of the oceans, mostly along the continental shelf (6). 

In shallow water the SOFAR channel rarely exists, and 
the effects of both bottom and surface reflections are signi- 
ficant. Some typical linearized velocity profiles are shown 
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in. figure 5* f’or uniform velocity gr-adients the ray path 

■the. 

equations as described for^surface channel may be used for 
the bottom reflection channel path calculations. For a, b, 
and c of figure 5» these calculations vjould have to be made 
in several steps. 

Except for the addition of a bottom reflection loss 
the bottom reflection channel losses are the same as those of 
the surface channel. Because of the greater channel depth it 
is obvious that the spherical spreading loss will extend for a 
greater distance. ■‘■his will be computed as R vice R/2 as^the 
surface channel. In general, experimental data shows that 
the total losses for this type of channel approximate spherical 
losses for the whole transmission distance, •‘•’he surface channel 
has total losses somewhat lov^er than spherical spreading. 

The bottom reflection loss is a complex function of the bottom 
material and its structure. For grazing angles less than 
fifteen degrees hov;ever, the losses are small. A good expe- 
rimental figure is 0.3 db per bounce (7). For larger grazing 
angles the loss becomes quite large. Figure 6 is a transmis- 
sion path diagram of the bottom reflection channel. Note that 
rays which graze the bottom at about I 5 degrees, leave the 
source at a much smaller angle. 

Sample calculations are made for both the surface 
channel and the bottom reflection channel, in Appendix A. 

Several computations are plotted in figure 7 together with 
some experimental data for comparison. In all cases the 
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Figure 5. Shallow water velocity profiles. 
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Figure 6 



The bottom reflection channel 
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Figure 7. Sound pressure level at 50 miles as a function 
of frequency. All data corrected to a source level 
• of 179 db. at a radius of one meter. References for 
experimental data in parentheses . 
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sources vrere corrected to a sound pressure level of I 85 db 
(ref. 0.00002 newt./m 2 ) at a 1/4 meter radius. 

SPL = 20 log^^Q P 

Po 

This is low enough to permit simultaneous transmission over 
each of ten channels. 

4. Noise . Ambient sea noise is a complex array of processes 
which are not all fully understood (13). Most noise sources 
are very frequency dependent. Belovr 1000 cps ship and traffic 
noise is the most significant. It usually peaks between '^00 
and 800 cps and falls off rapidly at higher frequencies. For 
all frequencies of interest above 1000 cps the ambient noise 
is primarily dependent on sea state for all but calm x-reather. 

The primary sources of this noise are thought to be wind, spray, 
cavitation, and precipitation. Other sources are marine life, 
molecular and thermal agitation and earthquakes. Table 1 
gives the average relations between wind force, vjind velocity, 
wave height, and sea state, figure 8 is a plot of the average 
ambient noise in the ocean versus sea state. The highest ex- 
pected values are about 3 above the average. A possible 
explanation is an irregular gustlng wind. It should also be 
noted that deep water noise levels are generally 5 lower 
than shallow x-iater levels. 

If a transducer is to be carried on a ship, the ship’s 
self noise is very important (14). At low speeds this will 
correspond to the ambient noise at sea state 2 , but it may 
contain some high level discrete frequencies. At speeds above 

/ 
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Figure 8, Ambient 
shallow water 



noise in the ocean(13). Specific 
data Indicated as follows: 
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(15) 

( 16 ) 
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about fifteen knots, cavitation noise can become much greater 
than the ambient noise. 



Table 1 



Sea 

State 


Beaufort 

wind 

force 


Wind Speed 
knots 

(raeters/sec . ) 


VJave Height 
feet 
(meters ) 


Sea 

Description 


2 


3 


7-10 

( 3 . 5 - 5 . 0 ) 


1-2 

( . 3 - .6) 


Scattered whitecaps 


3 


4 


11-16 

( 5 . 5 - 8.0) 


2-6 

( .6-1.8) 


Frequent whitecaps 


4 


5 


17-21 

( 8 . 5 - 11 . 0 ) 


5-10 

(1. 5-3.0) 


Many whitecaps 


5 


6 


22-27 

(11.5-14.0) 


8-17 

(2. 5 - 5 . 0 ) 


Whitecaps everywhere 


6 


7 


28-33 

(14.5-17.0) 


12-26 

(3. 5-8.0) 


Heaped up sea, 
blown spray 



5 . Source Power . The output power of a transducer is limited 
by cavitation and acoustic non-linearity . Cavitation occurs 
when the local pressure at or near the transducer is less 
than the vapor pressure of the water. A sample calculation 
in Appendix B shows that this limits the transducer to a 
sound pressure level of 195 (ref. 0.00002 newt./m2) at a 
transducer depth of twenty feet. 

The source level of the transducer is limited by the 
approximations required to obtain the acoustic wave equations 
from the equations of fluid dynamics. The governing ap- 
proximation comes from the following total derivative (1?): 

du _(ki . u<^u 
dt <)t <)r 
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u = fluid velocity vector 
t = time 

r = fluid particle position vector 
If the second term on the right can be neglected, then the 
wave equation can be satisfied. This is v/orked out in Ap- 
pendix B for spherical waves. The result is that the sound 
pressure level at one meter radius must be less than 254 db 
(ref. 0.00002 newt./m^), 

5. Random Process Channel i-iodel . An input-output relation 
for the hydroacoustic channel will now be developed. The 
parameters which will be considered are: signal energy, noise 

energy, frequency dispersion, time dispersion, and multipath 
reception. 

Since the ambient sea noise is due to a large number 
of small sources, it will be assumed to be a zero mean, 
Gaussian random process with spectral density N^. Nf is then 
the noise energy density, -^his noise is independent of the 
message sent. Therefore the received signal will be a 
message filtered by the channel plus additive Gaussian noise. 

This channel filter function is unknown however. To 
arrive at a channel description, tvro assumptions will be 
made. First, along any given transmission path, the channel 
is a pure attenuator. Second, the receiver receives the 
signal over a large number of paths, with the received phase 
for each path being a uniformly distributed random variable. 
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Wozencraft and Jacobs show that this description leads to 
a random process channel model, v;ith the following input- 
output relationships (18): 

s (t ) = m(t ) cos input 

r(t) = a(t)m(t) V? cos (Wot-&)+n(t ) putput 

in which 

p(a) = 2ae“^^ a— 0 

p(4) = 1 0£O‘^2tt 

2tt 

The amplitude a(t) is thus a Rayleigh-distributed random 
variable. Experimental data verifies this model in many 
instances, but other distributions have also been observed 
( 19 * 20 ), For the work to folloxi, a(t) x^rill be assumed to be 
Rayleigh -distributed. 

Before a communications system can be developed, 
frequency and time dispersion must be discussed. i'Vequency 
dispersion will limit the baud length.. Experimental data 
shows that in general, for ranges of interest, frequency 
dispersion is less than 2 cps (21). It is expected that the 
system will use a ba\id corresponding to a much larger band- 
width than this. 

Time dispersion of a signal is the range of delays 
observed over paths of different length. It limits the rate 
at which signals may be sent over conventional systems. 

■^'here is very little experimental data on shallow water time 
dispersion. The calculations in Appendix C show that about 
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1/5 second is the longest time that need to be considered, 
at a range of fifty miles. This is about 40 percent of the 
dispersion for the SOFAR channel I'lhich is severely limited 
in this respect (22), and it corresponds to an experimental 
value observed in the Florida Straits ( 23 ). 
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Part II 

COKIiUNICATIONS SYSTEM 

1. Frequency Selection , The optimum frequency is the one 
which allows the greatest bandwidth and the signal to noise 
ratio. Neglecting the noise, low frequencies are the best. 
Spreading losses are about the same for all frequencies. 

For a state 3 sea, surface reflection losses are 0.5 db./ 
bounce at 800 cps, 0.95 db. /bounce at 1000 cps, and 2.0 db./ 
bounce at 1500 cps. This becomes ;iorse for higher frequencies 
or higher seas. Attenuation also increases with frequency 
and is first significant for the fifty mile range at about 
2000 cps. Therefore, primarily because of surface losses, 
an upper limit of 1,5 kc will be used. A lov/er limit of 800 
cps will be used. At low frequencies it is possible to get 
strong harmonics from shipboard machinery. T'^is is general- 
ly worst at about 500 cps. One prime example of low frequency 
noise is propeller harmonics vjhich may be very strong, 

A plot of noise and fifty mile signal strength is 
shown in figure 9 as a function of frequency. 'I’his plot 
uses the same data as figures 7 and 8. A noise level cor- 
rected to a 10 cps bandwidth is also shown, ^his is suf- 
ficient for the correlation time of 1/6 second which will 
tentatively be used. Experiments show that at^ong ranges, 
the signal will remain about 80^ coherent by autocorrelation 
for up to 1/4 second (23 )• 
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SPL -(dk) 




Figure 9. Expected signal-noise separation in a state 
3 sea, (See figures 7 and 8.) 



/ ^ 
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2* Binary Sl.frnallln.g- , A simple binary communications 
system will now be developed for the shallow water channel. 
The signal set will be an orthogonal pair of cosine vxaves , 
as f ollovfs ; 

So(t) = V? cos A'qP 

Si(t) = -V? cos A>^t 

This set is considered the easiest to use for this system. 
Signals which are time orthogonal are not well suited to this 

channel since the time of arrival is unknown. Since both 

phase and amplitude are random processes, neither phase nor 
amplitude modulation was considered suitable. 



Vihen the signal So(t) is sent, the received signal is 
of the form: 

r(t) = ^ ai cos(^^ot+©i) + n(t) 

1=1 



where the i's represent the different multipaths. As noted 
in Part I, this leads to the following random process model: 
r(t) = a(t) cos (W o'*' + 9) + n(t) 
in which a and 9 are random variables with the following 
probability density functions at any instant of time: 
p(a) = 2a e”^^ 0 ^ a 

p(0) = 1/27T 0 ^ 9 ^ 27T 

This model and the follovring assumptions will be used to 
arrive at a suitable binary receiver: 

a. a(t) and 9(t) are constant over the signal inter- 
val T, but vary randomly from interval to interval. 

b. n(t) is a white Gaussian noise random process, 
of spectral density No, and it is independent of the signal. 
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c. The two signals, So(t) and S 2 (t) have equal 
energy and are equiprobable , 

d. There is no intersyrabol interference. 

A maximum likelihood receiver will be used. For this 
receiver, that message is selected for vrhich the probability 
of the received signal conditioned by the messages is highest. 
The details of this receiver have been vrorked out by Turin 
( 26 ), and for more general time-varying cases by Kailath (2?) 
and Price (28). In the case of the narrow band signals chosen 
for the signal set, the receiver was determined to be a de- 
vice which correlates, squares, and sums the quadrature com- 
ponents before making a decision. The correlation receiver 
is shown schematically in figure 10. Figure 11 is an envelope 
detector receiver which can be shown to be essentially the 
same . 



The correlation receiver can be explained heuristic- 
ally by a geometric argument. The quadrature noise com- 
ponents of the received waveform are zero mean Gaussian ran- 
dom variables of variance Nq/ 2, where Nq is the mean square 
noise power. The quadrature signal components of the received 
waveform are also zero mean Gaussian random variables, but 
with a variance of Eq/ 2, where Eq is the mean square signal 
power. If a signal is present the received quadrature com- 
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Figure 10. 



Binary correlation receiver. 




Figure 11. Binary envelope detector receiver. 
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ponents are therefore zero mean Gaussian random variables 

with variance N , . These Gaussian random variables are 

2 

illustrated in figure 12, in the So -plane and the s, -plane 
for the case when So(t) is sent. In each case the amplitude 
of the resultant vector is a Rayleigh distributed random 
variable . 

The error probabilities for the likelihood receiver are 

developed in Wozencraft and Jacobs (29), by computing it 

first for the random phase case only, to obtain: 

-R/2N 

P (£) = i: e 

Averaging this error over a Rayleigh amplitude distribution, 
the following error probability is obtained for the random 
phase and random amplitude channel: 

~ 2 + E/N 

Using this formula, the efficiency of the simple binary 
channel may now be determined. Referring to figure 9, v/ith 
a 10 cps bandv^idth and the weak signal line, table 2 summarizes 
the binary channel error probability at various frequencies. 

Table 2 



Frequency 


E/N (db. ) 


e/n 


P{£) 


800 


29 


790 


0.0013 


1000 


23 


200 


0.0050 


1250 


15 


51.5 


0.030 


1500 


8 


6.3 


0.12 



3. Channel Capacity and Coding . In 19^3, C.E. Shannon 
showed that whenever signals are transmitted at some rate, 
R, less than channel capacity, C, then an arbitrarily small 
error rate can be achieved by coding with a sufficiently 
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Figure 12. Received signal In the -plane and s, -plane 
when message s^Ct) is sent. 




Figure 12A. The binary symmetric channel. 
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lon^ block length ( 30 ), For the binary symmetric channel 
with crossover probability p, as shovjn in f igure 12A, the 
channel capacity has the rather simple form ( 31 ): 



where T is the signal time required to achieve p. It can 
also be shown that for n arbitrary channels being used in 
parallel, the capacity is: 



-the 

where/^Ci are the individual channel capacities. This last 
formula holds for simple parallel use in vJhlch a decision 
is made on each channel. If the decision is made on the col- 
lective output probabilities of the individual channels, then 
the channel capacity is somewhat greater, since in the first 
method some information is lost by the intermediate decision, 

A sixty word per minute data rate requires approxima- 
tely thirty bits of information per second. If an average 
channel capacity of five bits per second can be obtained, then 
a minimum of six channels are required for effective trans- 
mission. At transmission rates very close to capacity how- 
ever, very elaborate coding schemes are required to achieve 
a reliable transmission. In general, it can be shown that 
for proper coding the probability of error is bounded by 
functions of the transmission rate ( 32 ): 




bits per second 



n 




e-n(2L(R) + 0(n) ^ p^^j^^-nE(R) 



where E^^R) and E(R) depend on the channel, and n is the code 
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block length. For rates less than channel capacity, both 
E^CR) and E(R) are positive. These functions will not be 
evaluated for this channel. 



For this fading channel, the rate is limited by the 
number of parallel channels. The signal interval cannot be 
shortened to much less than one sixth of a second because of 
intersymbol interference. The number of channels is limited 
by the useable frequency band. Above about I5OO cps, very 
little cfipacity is gained for the addition of an extra channel. 
Assuming a signal interval of one sixth of a second, and trans- 
mitting five information bits during each interval, the ef- 
fect on the rate as a function of capacity of adding more 
channels is summarized in table 3 » and figure I3 . A minimum 
useable frequency of 900 cps and a frequency separation of 
20 cps was assumed. 







Table 3 "' 








Channel # 


frequencies 


P(^) 


AC 


C 


R/C 


L 


900 & 920 


.0028 


5.8 






2 


940 & 960 


.0040 


5.8 






3 


980 & 1000 


.0057 


5.7 






4 


1020 & 1040 


.0076 


5.6 






5 


1060 &1080 


.0100 


5.5 


28.4 


1.06 


6 


1100 & 1120 


.0124 


5.4 


33.8 


.89 


7 


1140 & 1160 


.0154 


5.3 


39.1 


.77 


8 


1180 & 1200 


.0213 


5.1 


44.2 


.68 


9 


1220 & 1240 


.0263 


4.9 


49.1 


.61 


10 


1260 & 1280 


.0333 


4.7 


53.8 


.56 


11 


1300 & 1320 


.0408 


4.5 


58.3 


.51 


12 


1340 & 1360 


.0562 


4.1 


62.4 


.48 


13 


1380 & 1400 


.0685 


3.8 


66.2 


.45 


14 


1420 & 1440 


.0918 


3.3 


69.5 


.43 


15 


1460 & 1480 


.114 


2.9 


72.4 


.41 



* Based on a power level allovjlng ten channels. For less 
channels is better, for more worse (See Part I). 
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Figure 13. The effect of adding additional higher frequency 
channels on capacity and rate. (Transmission at 30 
bits per second.) 
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There are several Interesting methods of coding and 
decoding which will be discussed after developing a com- 
munications system for the hydroacoustic channel. In the 
next section, the system will be developed first using simple 
transmission and then diversity and block coding. Block 
coding is taking a block of input digits and coding them into 
a longer group of coded digits for transmission. The error 
will be calculated for a decoder which makes an intermediate 
decision on each bit. As noted this is not the most efficient 
use of the channel, but it is considerably easier for decoding 
and for error calculations. 

4. Representative Communications Systems . A simple system 
vfhich can transmit the necessary 30 bits per second is com- 
posed of five parallel binary channels. Each one sixth of a 
second, a five bit symbol is transmitted, one bit per channel. 
To avoid shipping noise a minimum frequency of 900 cps x^lll 
be used. Table 4 lists the individual crossover probabili- 
ties for these five channels. The orthogonal signal pair are 
12 cps apart to allow a 1/6 second signal interval. The 
channels are separated by 20 cps, to permit easier tuning, 
in case of a doppler shift which could be as high as 15 cps 
for two ships at 15 knots. 



Channel # 


Table 4 
freauencies 


P(^) 


1 


900 & 912 


0.0014 


2 


932 & 944 


0.0019 


3 


964 & 976 


0.0024 


4 


996 & 1008 


0.0030 


5 


1028 & 1040 


0.0038 



s. 
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The error rate per five bit symbol for this transmission 
system can easily be computed as: 

P(<f) = 1 - P(no bit errors) 

For independent bit errors the error probability is easily 
computed to be O.OI 3 per five bit symbol. Due to receiver 
imperfections, and variation of the channel from the model 
vfhlch will be discussed in the next section, this can only 
be considered an order of magnitude. Therefore, methods of 
improving the reliability will be investigated. Shannon’s 
coding theorem discussed in the last section indicates that, 
the error rate may be improved without an increase of power. 

The easiest method to implement for this channel is 
not normally thought of as coding, though it is in reality 
a simple code. This is frequency diversity. For illustration 
fifteen parallel channels will be used, with each bit being 
transmitted over three channels instead of one. T© maintain 
the same transmitter power, each channel will use only one 
third of the power as for the simple system binary channels. 
The error probabilities for the fifteen channels are listed 
in table 5* The frequencies are spaced as for the simple 
five channel system. 



Table 5 



Channel # 


7{d) 


Channel # 


P(^) 


T 


.0025 


0 

y 


.0264 


2 


.0059 


10 


.0324 


3 


.0074 


11 


.0385 


4 


.0090 


12 


.0446 


5 


.0117 


13 


.0550 


6 


.0149 


14 


.0705 


7 


.0184 


15 


.0840 


8 


.0224 
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A block diagrani of the triple diversity receiver is shov/n in 

figure 14. The are weighting functions which depend on 

the signal-noise ratio for the particular binary channel (33): 

Wj_ = Bj/Nl 

l+3j_/Ni 

For the fifteen channels being used, they rs-nge from 0.99 
for channel number one to 0.86 for channel number fifteen. 

For a simple equal ene?rgy L-fold diversity system the error 
probability can be computed from the following formula (3^): 




This is plotted for L equal to 1, 2 , and 3 i'd figure 15. As 
an example, if p is 0.04, three-fold diversity yields a bit 
error probability of O.OOO 65 . For a five bit symbol this 
vjould become an error rate of 0.0033. Assuming that p equal 
to .04 is a good average from table 5i this is a four times 
improvement over simple five channel transmission. A general 
error calculation method for unequal energy and noise on the 
different diversity channels is worked out in Appendix D. 

If data is to be transmitted X'lith an accuracy better 
than this diversity system, coding will be required. One 
method, block to block coding will be discussed now. Other 
methods vrill be compared to it in the next section. A prac- 
tical block to block coding system is the parity check code. 
For illustration, each five bit symbol will be encoded into a 
fifteen bit code word. The first five bits of the code vjord 
are the five information bits; the remainder are modulo two 



A 



31 



r (t) 







Figure 14 



Three-fold diversity receiver 
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Figure 15. Binary receiver error probabilities for 
the Rayleigh fading channel. 
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sums of several of these five bits. Proper selection of the 
code, and use of an appropriate receiver, a-ifill permit cor- 
rection of any three bit or lesser error in the reception 
of the fifteen bit code word (35 )• For this code the error 
probability then becomes ; 

P() = l-p(no errors) - p(one error) - p(two errors) - p (three 

errors ) 

Carrying out this calculation for the data of table 5» the 
error probability is approximately 0.004 per five bit symbol, 
or approximately the same as for triple diversity. Further 
improvement can be accomplished by coding by use of a longer 
block length, for example, coding tvio five bit symbols into 
a thirty bit code word, or three symbols into a 45 bit vjord. 

5. Suitability of Coding- Jochnioues . In the last section 
simple transmission, frequency diversity and a short block 
length parity check code were compared, P’nere was a signi- 
ficantly smaller error using diversity or the coding scheme. 

To achieve a still smaller error probability, much larger 
block length or other types of codes would be required. 

Various coding techniques I'fill now be discussed, xfuich may 
be used if the added complexity can be justified by the im- 
proved performance, 

A block to block parity check code was used as an 
example in the last section. For this type of coding system, 
the number of computations required by the decoder increases 
roughly as the cube of the block length (3^), and as noted. 
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only specific error combinations can be corrected. This 
limits its capability of correcting a large burst of errors, 
except for very long codes. In general, block code receivers 
make a decision on each bit before decoding. This means that 
the confidence level in each bit is lost. The perf ormance 
can be Improved by finer quantization in the receiver before 
decoding. It can be shown, however, that three level quanti- 
zation will achieve half of the possible improvement, and 
that finer quantization than this gives an increasing small 
improvement (37 )• Three level quantization is best described 
by the binary erasure channel shown in figure l6. For this 
channel the bias level must be properly set so that both p and 
q are small, where p is the crossover probability, and q is 
the erasure probability. This type of channel is more dif- 
ficult to instrument than the simple binary channel and has 
not been used as extensively ( 38 ). 

Some of the dravjbacks of block to block coding may 
be met with convolutional codes (39) • A simple encoder is 
shown in figure 17 (40). Each transmission interval, v 
parity checks on a k bit register are transmitted. Between 
intervals, one bit is moved out of the k bit register on the 
right, and one nevj information bit is moved in at the left. 
Each information bit is in the transmission sequence for k 
intervals, and there are an average of v transmitted bits per 
infolmation bit. Since at least five information bits must 
be transmitted on fifteen channels six times per second, v is 
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Figure 



iniorrna.'iien bUs 




16 . The binary erasure channel. 
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Figure 1? 



Convolutional encoder 
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limited to three for the channel being described. 

Two good decoding schemes are known for convolutional 
codes: sequential decoding (4l), and threshold decoding (42), 
In the former an n bit message is treated as a tree with n 
sequential nodes. A very good reliability can be achieved 
if a computer with a large enough memory to backtrack through 
many paths along the tree is used. This may be 2^^ bits or 
more (39). Once an error is made because backtracking is no 
longer possible, a large block of errors may be made. This 
is satisfactory in a feedback system where repeats may be 
requested. Since the transmission time over a fifty mile 
hydroacoustic channel is approximately a minute, this type of 
channel is best suited for broadcast or prearranged time 
sharing, vrith little or no possible feedback. A simpler 
decoding scheme for convolutional codes which avoids these 
block errors at the expense of a higher average error rate 
is threshold decoding. 

The threshold decoder uses only shift registers, modulo 
two adders, and a threshold device. For good results each 
bit must remain in the register for more than one fifteen 
bit transmission time, since for this short coding length 
the results are about the same as for triple diversity as 
noted on figure 15 (43). For short block lengths, the results 
viith threshold decoding on a convolutional code are somewhat 
better than for block to block coding, vrith no more complexity. 
This is because it can take advantage of the a posteriori 
probabilities in the binary channels. 
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Part III 
DISCUSSION 

1. General . In this section the assumptions leading to 
the channel model and the resulting communications system 
v:ill be discussed first. The Importance of the assumptions 
to the resulting comm.uni cat ions system, and their expected 
variations vjill be brought out. Next the hydroacoustic con- 
ditions in various ocean areas at different times of the year 
will be discussed. 

2. xA.s sumptions . The most important assumptions made in the 
preceding development have been the follov;ing; 

a. There is no intersymbol interference. 

b. The multipath reception produces a Rayleigh 
amplitude distribution. 

c. During each signal interval, the received phase 
and amplitude are constant. 

d. The receiver is synchronized to the transmission 
interval. 

e. The noise is vihite and Gaussian. 

f. Errors on different binary channels are independent. 

Assumptions a., c., and d. are related. To minimize inter- 
symbol interference, the signalling interval must be long 
enough for most of the energy transmitted during a given 
interval to be received during an equal delayed interval at 
the receiver. The time used is based on the calculation in 
appendix G. The intersymbol interference can be treated as 
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an additional noise component which increases the crossover 
probability, or if it can be estimated, used in a receiver 
with memory. By making the signalling Interval long enough 
for intersymbol interference to be negligible, the interval 
may become too long to assume a constant phase and amplitude 
during the interval. This problem can be corrected by using 
time diversity reception during each signal interval. Since 
a small amount of intersymbol interference is unavoidable, 
and also because the transmission time betvreen the terminals 
will vary, it will be impossible to obtain perfect synchroni- 
zation at the receiver. Range changes may be corrected by 
a doppler circuit. Transmission path variations should be 
fairly small so that the time variations due to this phenomena 
vii 1 be small compared to the signal Interval. In this case 
these variations are also correctable. In summary, these 
assumptions can be met as vrell as desired by proper selection 
of the signalling interval and the use of time diversity 
reception if required. 

Experimental data indicates that the Rayleigh ampli- 
tude distribution is a good assumption much of the time, but 
that it may differ (19,20). For long ranges these differences 
are most evident in the higher moments, and do not differ 
significantly in the mean or mean square. At short ranges 
only a few paths may be of interest, and in this case the 
distribution is more favorable than a Rayleigh distribution 
(19). Figure 18 compares a typical short range amplitude 
distribution vrith the Rayleigh distribution. 
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Figure 18. Typical amplitude distributions. 




Figure 19. Effect of a sloping bottom on communications 
via the bottom reflection channel in shallow water. 
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Although most ambient sea noise can be considered 
x-Jhite and Gaussian for the narrox'i band binary channels, other 
forms of noise may be present. Rotating machinery harmonics, 

6r marine life may affect only one or txMo frequencies. Kost 
of these noises are most severe under 1000 cps, but a source 
on a receiving ship is potentially dangerous. A careful noise 
analysis of potential receiving vehicles will permit optimum 
frequency selection with respect to the noises. Impulse 
noises may also be present. These can be minimized by clip- 
ping the signal at a level which depends on the average im- 
pulse noise pov^er and the expected signal levels. 

In arriving at ;ji system error probabilities, errors 
were considered independent of other channels. If these 
binary errors are not independent, then the resulting system 
error probabilities x^ill be different. If for example the 
errors on tvjo of three diversity receivers are correlated, 
the resulting bit error is only about as good as double and 
not triple diversity. No good data is available on how closely 
correlated transmission over adjacent frequencies is. 

3 . Ocean Conditions and the Channel Model . Several factors 
affect the prevailing ocean hydroacoustic conditions (40). 

The temperature, and hence the vertical profile of the velo- 
city of sound vary with both latitude and time of year, and 
in some cases the time of the day. The depth of the water 
and the bottom profile affect bottom reflection transmission. 
Vfeather conditions affect both noise and surface reflection 
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losses. All the calculations made have been for a moderate, 
state three, sea. Heavier weather conditions reduce the 
reliability considerably. 

Surface trapped channels occur only when the water 
temperature is isothermal or perhaps increasing vxith depth. 
4Viis condition is somewhat unstable and can occur only when 
the surface vraters are being cooled. Experimental data 
indicates that these conditions exist most of the year in 
polar regions, during winter in mid-latitudes, and only 
rarely in tropical waters. In mid-latitudes the layer depth 
will be from 1000 feet to I 50 O feet in the winter, but a 
sharp thermocline will exist near the surface in the summer. 
Since water depth is not important to the surface channel, 
the model developed will apply to all high latitude ocean 
areas in the winter, but only to shallovf water, or polar 
waters during the summer months . 

The bottom reflection channel can also be quite 
variable as to results. For any given depth, the bottom 
s].ope is quite important, as is illustrated in figure 19. 
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Part IV 

CONCLUSIONS AND RSCONDiENDATIONS 

1. Conclusions « In Part I it was shown that the signal- 
noise ratio for the surface channel and the shallow water 
bottom reflection channel could be estimated. It was further 
shown that the appropriate model was a Gaussian random pro- 
cess, in which the received signal had a Rayleigh amplitude 
distribution and unknown phase with additive vrhite Gaussian 
noise. Using this model an appropriate binary receiver was 
found to be a quadrature correlation device. A combination 
of five of these receivers could be used to send data at 60 
wpm vrith a moderate error rate. Improvement in the error 
rate vias best obtained by frequency diversity, and could be 
further improved by coding at the expense of complexity and 
cost. In Part III, it was shovm that the model v:as reasonable 
for shallow water any time if the bottom did not slope away 
from the transmitter too fast, for mid-latitudes in the 
winter, and for polar regions most of the year. 'P’nese 
results are applicable for moderate weather conditions, and 
ranges of up to fifty miles. 

2. Recommendations . Although the basic formi of the receiver 
could be developed, the actual parameters could only be 
estimated. A better understanding of these parameters is 
needed for the actual development of a suitable system. 

Some items of Interest which warrant additional experimental 
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research are the following: 

a. the noise spectrum of potential receiving 
vehicles . 

b. the range and distribution of multipath time 
delays for various channels. 

c. the length of time during which phase and ampli 
tude of a multipath signal I'Jill remain coherent 

d. the correlation of errors on binary channels of 
adjacent frequency. 



Ar.-EMDIX A 



1 . Surface Channel Sample Calculations . 

For these calculations the following assumptions will be made 
Isothermal water with a velocity gradient of O.CO 52 m./ 
sec . /m. 

Frequence ... 1 kc. 

Hange ... 50 miles. 

layer depth ... 30*5 meters. 

Sea state ... 3 

Source ... I 85 db (ref. 0.00002 nev:t./m,) at l/'4 m. 
radius. This is 1?3 db at a 1 meter radius. 

a. H = C2 ^ 

a c os G, 

0=0 

c^ = 1500 m/sec . 
a = 0.0052 m/sec/m. 

Y = 30.5 m 

. . - a 

COS 0„ = c, 1-Q« 

c, +aY 2 

0 « = 0.0152 radians = O .87 degrees 
R =4380 meters 

b. Spherical spreading loss. 

L = 20 log^^R/2 = 20 log^^2190 = 6 ? db. 

c. Cylindrical spreading loss. 

L = 10 log 50 miles = 10 log 36.7 = 15*7 db. 

R/2 

d. Surface reflection loss. 

Number reflections = 50 miles = 9«2 

2 R 
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Number reflections = 10 

Loss per bounce = 0.95 <3.4 for a 4 foot wave height 
L = 0.95 X 10 = 9.5 <ib. 

e. ^i-bsorption. Negligible for this range and frequency. 

f. Total losses. 6 ? + 15*? + 9*5 = 92 db. 

g. Signal level. 1?3 - 92 = 8 l db (ref. 0.00002 nev:t./m.) 



2 . Bottom Reflection Channel Calculations . 

The same assumptions are made as for the surface channel 
except there is no layer, and; 

VJater depth ... 300 meters, flat bottom. 

Velocity gradient ... 0.0?5 m. /sec/m. 



a . R = 



a cos Oo 



(sin ©,) 



R = 3500 meters 

b. Spherical loss. 

L = 20 log^^3500 = 71 db. 

c. 'Cylindrical loss. 

L = 10 lov 50 miles = 10 log,„23 = I 3.6 db. 

3500 

d. Bottom Loss. 

Number reflections = 50 miles = 11.5 

2R 

Number reflections = 12 

L = 12 X 0.3 = 3.6 db. 

e. Surface loss. 

L = 12 X 0.95 = 11.^ (i4. 

f. Total loss. 71 + 13.6 + 3.6 + 11.4 = 100 db. 

g. Signal level . 173-100 = 73 <i4. (ref. 0.00002 newt./m.) 
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1 . Cavitation . 

Cavitation occurs v;hen the acoustic pressure plus the ambient 
pressure is less than the vapor pressure of water. The 
following nunerical values will be used in this calculation: 

vapor pressure = 0.6 psi at 85 degrees P. (decreases vrith T) 
atmospheric pressure ... 14.7 psi 
water depth =20 feet 

a. Static pressure. 

p = 14.7 + 0.444 X 20 = 23.5 psi 

b. Allov;able acoustic pressure. 

P = 23.5 - 0*6 = 22.9 psi = 1.58 X 10 nevit./n, 

c. Sound pressure level. 



Allowing 3 db for peak versus mean square, the S?L is 
limited to 195 db (ref. 0.00002 nevjt./m^) 

2 . Non-Linearity . 

For the acoustic approximation v.'e require; 

uJu ^u 
i r ^ J t 

Assuming spherical harmonic waves, the following equations 
apply; 

P = e< 5 (*'*t-kr) 



SPL =20 log,^p/p(ref ) 
p(ref) = 0.00002 newt./m^ 

SPL = 20 log ,/?. 9 X 10’ = 193 db 



r 
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Taking the indicated derivatives we nov: obtain; 

v> [l. + jk 1 « 0 "* 

jd [r^ r -I 

let f, = = 1000 cps 

r = 1 meter 

then, since k = , the bracketed term on the left is less 

than 10. '•‘•‘he Inequality can be simplified to 

£• « 

e. 

In sea X'Jater, 1000 kg/m. 

2 , i ^ 

o). = 4 TT" }C 10 at f, = 1 kc. 

so that the inequality becomes: 

P « 4 X lo’ 

8 

Let: ting p = 10 , vie obtain the folloviing sound pressure level 
SPL = 20 log^^p/p(ref ) 

= 20 log 10*" 254 db (ref. 0.00002 nevrt./m^. ) 

'*0.00002 
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Ax-.-EJIDIX C 

1. •‘•'iFe i)lsT:>ers ion . The Kaximum difference in the length 

of time it takes a signal to reach the receiver via two dif- 
ferent multipaths limits the rate at which signals may be sent, 
‘^his difference is due to sound velocity fluctuations, and 
path length fluctuations. This time difference vrill now be 
estiraated for the special case of a bottom reflecting channel. 
Variations xfill be less severe for surface channels or for 
shallovier water. It will be assumed that the average velocity 
of sound differs by a maximum of four meters per second over 
different paths. This corresponds to an average temperature 
difference of 1.0°P. This is a reasonable value since va- 
riations over several miles are normally quite small in open 
vjaters. It will further be assumed that the ray paths are 
straight lines. This is illustrated in figure 20. Since 0 
is limited to 10° - 15 ° at the bottom for good reflection 
coefficients, an average value of 6° will be chosen. This 
is reasonable since the most significant paths are those 
leaving the transducer X’Jlthln a fev: degrees of the horizontal. 
All these paths will reach the bottom at approximately the 
samie angle, so that the average variation viill be only about 1°. 

Using these values, a range of 50 miles, and a sound 
velocity of I 5 OO m/sec., the time variation will now be 
computed . 

T = D/V = slant range/velocity 

R = D cos 0 = horizontal range 
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clepih 



Figure 20. Model for time dispersion calculation 



50 



” ~ dv - H tan QdO 

V* V 

B dv - R OdS 

V 

where 

H = 50 miles = 80,500 meters 
V = 1500 ra/sec. 
dV = 4m/sec. 

0 = ,104 radians 
dO = .0174 radians 

Carrying out the calculation the time variation is 0,240 sec. 
This is a variation from the straight line path, however. 

Sven the shortest paths will range from the straight line by 
.04 - .06 sec. This leaves a total variation of about 0.2 
seconds vjhich compares favorably with the only experimental 
data available (23). 



: ^ 
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ArPSNDIX D 



1. Diversity Error Calculation . In part II it is noted 
tdat the diversity receiver makes a decision that iCo is 



sent if: 



2L ^ 2L ^ 

Where the Slj. are statistically independent zero mean Gaussian 
random variables. Since the vij. are all approximately equal 
they will be neglected and the error probability vjill be 
computed for making the decision on: 




By pairwise addition of the quadrature components on each 
channel, the decision can also be represented as: 




vrhere the bj and the b^ are statistically independent Rayleigh 
distributed random variables, x>rith the follovring probability 
distributions . 



P(bj 



4bj 

Ej +Nj 



exp 




E . 



0 < hj < ^ 



P(b; ) = ^ exp - 2bi. 



0 be < 



The error probability can be easily calculated if the pro- 
babilities of the two sums are known. These will now be 
■out in a closed form. Define a nev? variable: 



kj = bj 

then it is easily shown that 



p(k. ) = 2 

TTX 






_ 2k; 

E . -I-N,. 



J 



exp 
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The probability distribution of the sum of the kj can be 
found by noting that the characteristic function of the sum 
of independent random variables is the product of their in- 
dividual characteristic functions. The characteristic fun- 
ction of k is : 



fpikj ) e dkj = 2 (_1_ \ 

2; -rlvt s+2/Ej +N/ 

Therefore the characteristic function of the sum of the kj is: 

[p(^kj= R^.) e-sRj dR: = TT ( — 2 — j (__1_ 

ji/ Ej +Nj ^ s+2/Sj ^ 

By a partial fraction expansion it is easily shovrn that the 
probability density function of the sum of the kj is the 
sum of L exponentials . 



Knowing the probability functions, the probability of 
error can easily be calculated. If m<, is sent, and 



d=l 



k , = R 

do 



o 



a =1 



then the error probability is : 
p(^) = fp{R/ ) [ p(Ro) dRodR, 

-'o 

This integral can be solved in general, for the special case 
of tvjo-fold diversity it takes the folTowing form; 



p(f) = 1 - 



K.-K, 



K. 



1 



K, -Kj/ iK^-rX^ 



K, -kJ 



1 - 



K,+Kv/i 



53 



v:here : 






IC, = 

K3 = 

K, = 

j-'his last result is restricted to cases 'ut.ere K, T^ICjand 

If one or toth of these pairs were equal, the result 
would fee found in a similar v:ay, but would be slightly dif- 
ferent. Figure 21 shov/s the effect of varying signal-noise 
ratio on tvjo-fold divers itjr. 
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Figure 21. Diversity error with unequal signal-noise ratios. 
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